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Abstract. We complete the existing literature on the kinetic theory of systems with 
long-range interactions. Starting from the BBGKY hierarchy, or using projection 
operator technics or a quasilinear theory, a general kinetic equation can be derived 
when collective effects are neglected. This equation (which is not well-known) applies to 
possibly spatially inhomogeneous systems, which is specific to systems with long-range 
interactions. Interestingly, the structure of this kinetic equation bears a clear physical 
meaning in terms of generalized Kubo relations. Furthermore, this equation takes 
a very similar form for stellar systems and two-dimensional point vortices providing 
therefore a unified description of the kinetic theory of these systems. If we assume that 
the system is spatially homogeneous (or axisymmctric for point vortices) , this equation 
can be simplified and reduces to the Landau equation (or its counterpart for point 
vortices). Our formalism thus offers a simple derivation of Landau-type equations. We 
also use this general formalism to derive a kinetic equation, written in angle-action 
variables, describing spatially inhomogeneous systems with long-range interactions. 
This new derivation solves the shortcomings of our previous derivation [P.H. Chavanis, 
Physica A 377, 469 (2007)]. Finally, we consider a test particle approach and derive 
general expressions for the diffusion and friction (or drift) coefficients of a test particle 
evolving in a bath of field particles. We make contact with the expressions previously 
obtained in the literature. As an application of the kinetic theory, we argue that 
the relaxation time is shorter for inhomogeneous (or high-dimensional) systems than 
for homogeneous (or low-dimensional) systems because there arc potentially more 
resonances. We compare this prediction with existing numerical results for the HMF 
model and 2D point vortices. For the HMF model, we argue that the relaxation time 
scales as N for inhomogeneous distributions and as e N for permanently homogeneous 
distributions. Phase space structures can reduce the relaxation time by creating some 
inhomogencitics and resonances. Similar results are expected for 2D point vortices. 
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1. Introduction 

Kinetic theories of systems with long-range interactions are important to understand 
the dynamical evolution of the system and to study transport properties. They can 
be used to determine the timescale of "collisional" relaxation towards the Boltzmann 
distribution, in particular its dependence on the number of particles N. Furthermore, 
they can tell whether or not the system will truly relax towards statistical equilibrium. 
Indeed, it is not obvious that complex systems (like systems with long-range interactions) 
do eventually relax towards Boltzmann equilibrium because the ergodicity assumption 
which sustains the statistical theory can break down. Different kinetic theories have 
been developed over the years. 

Dilute gases: The first kinetic theory of many-particles system was elaborated by 
Boltzmann [1] for a dilute gas. In that case, the particles do not interact except during 
strong collisions. Boltzmann derived his famous equation from which he proved that 
the entropy increases (if -theorem) and that the system ultimately relaxes towards the 
Maxwell distribution of statistical equilibrium. 

Coulombian plasmas: Landau [2] derived a kinetic equation for a Coulombian 
plasma by starting from the Boltzmann equation and considering a weak deflection 
approximation. Indeed, for a Coulombian potential of interaction slowly decreasing with 
the distance as r _1 , weak collisions are the most frequent ones. Each encounter induces 
a small change in the velocity of a particle but the cumulated effect of these encounters 
leads to a macroscopic process of diffusion in velocity space. The treatment of Landau, 
which assumes that the particles follow linear trajectories with constant velocity in a 
first approximation, yields a logarithmic divergence of the diffusion coefficient for both 
small and large impact parameters but the equation can still be used successfully if 
appropriate cut-offs are introduced. A natural lower cut-off, which is called the Landau 
length, corresponds to the impact parameter leading to a deflection at 90°. On the 
other hand, in a neutral plasma, the potential is screened on a distance corresponding 
to the Debye length. Phenomenologically, the Debye length provides an upper cut-off. 
Later on, Lenard [3] and Balescu [I] developed a more precise kinetic theory that could 
take into account collective effects. This gives rise to the inclusion of the dielectric 
function |e(k, k ■ v)| 2 in the denominator of the potential of interaction appearing in 
the kinetic equation. Physically, this means that the particles are "dressed" by a 
polarization cloud. The original Landau equation, which ignores collective effects, is 
recovered from the Lenard-Balescu equation when |e(k, k • v)| 2 = 1. However, with 
this additional term, it is found that the logarithmic divergence at large scales is now 
removed and that the Debye length is indeed the natural upper lengthscale to consider. 
At about the same period, Hubbard [5] developed a test particle approach and derived 
a Fokker-Planck equation describing the relaxation of a test particle in a bath of field 
particles. He calculated the diffusion and friction coefficients by evaluating the first and 
second moments of the velocity deflection and took into account collective effects. It is 
interesting to note, for historical reasons, that Hubbard [5] was apparently not aware 
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of the works of Lenard [3] and Balescu [I] at the same period and that he developed 
his approach independently. However, the three approaches are closely related. Indeed, 
if we substitute the diffusion and friction coefficients found by Hubbard in the general 
form of the Fokker-Planck equation and perform minor transformations (a substitution 
that Hubbard has not explicitly made), one obtains the Lenard-Balescu equation! 

Stellar systems: In stellar dynamics, Chandrasekhar [TJ Ej developed a kinetic 
theory of stars in order to determine the timescale of collisional relaxation and the rate 
of escape of stars from globular clusters. To simplify the kinetic theory, he considered an 
infinite and homogeneous system. He started from the general Fokker-Planck equation 
and determined the diffusion coefficient and the friction force (second and first moments 
of the velocity increments) by considering the mean effect of a succession of two-body 
encounters^- Since his approach can take into account large deflections, there is no 
divergence at small impact parameters and the gravitational analogue of the Landau 
length appears naturally in the treatment of Chandrasekhar. However, his approach 
leads to a logarithmic divergence at large scales that is more difficult to remove in 
stellar dynamics than in plasma physics because of the absence of Debye shielding for 
the gravitational force. In a series of papers, Chandrasekhar & von Neumann pU] 
developed a completely stochastic formalism of gravitational fluctuations and showed 
that the fluctuations of the gravitational force are given by the Holtzmark distribution 
(a particular Levy law) in which the nearest neighbor plays a prevalent role. From 
these results, they argued that the logarithmic divergence has to be cut-off at the 
interparticle distance. However, since the interparticle distance is smaller than the 
Debye length, the same arguments should also apply in plasma physics, which is not 
the case. Therefore, the conclusions of Chandrasekhar & von Neumann are usually 
taken with circumspection. In particular, Cohen et al. [11] argue that the logarithmic 
divergence should be cut-off at the physical size R of the cluster (or at the Jeans length) 
since the Jeans length is the presumable analogue of the Debye length in the present 
context. These kinetic theories lead to a relaxation time of the form t^ ~ (N/hiN)t£, 
where to is the dynamical time. Chandrasekhar [12] also developed a Brownian theory 
of stellar dynamics and showed that, on a qualitative point of view, the results of 
kinetic theory can be understood very simply in that framework. In particular, he 
showed that the dynamical friction is necessary to reproduce the Maxwell-Boltzmann 
distribution at equilibrium and that the coefficients of diffusion and friction are related 
to each other by an Einstein relation. This relation is confirmed by his more precise 
kinetic theory. It is important to emphasize, however, that Chandrasekhar did not 
derive the kinetic equation for the evolution of the system as a whole. Indeed, he 
considered the Brownian motion of a test star in a fixed distribution of field stars (bath) 
and derived the corresponding Fokker-Planck equation^. This equation has been used 

| Later, Rosenbluth et al. |9j proposed a simplified derivation of the coefficients of diffusion and friction 
for stellar systems and for plasmas (without collective effects). 

§ Indeed, Chandrasekhar [7j [8] models the evolution of globular clusters by the Kramers equation 
which has a fixed temperature (canonical description) while a more relevant kinetic equation would be 
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to study the evaporation of stars from globular clusters [5J [T3J dH [J51 EE]. King [T7] 
noted that, if we were to describe the dynamical evolution of the cluster as a whole, 
the distribution of the field particles must evolve in time in a self-consistent manner so 
that the kinetic equation must be an integrodifferential equation. The kinetic equation 
obtained by King is equivalent to the Landau equation^. It is interesting to note, for 
historical reasons, that none of the previous authors seemed to be aware of the work of 
Landau [2] at that time0- There is, however, an important difference between stellar 
dynamics and plasma physics. Neutral plasmas are usually spatially homogeneous due 
to Debye shielding. By contrast, stellar systems are spatially inhomogeneous. The 
above-mentioned kinetic theories developed for an infinite homogeneous system can be 
applied to an inhomogeneous system only if we make a local approximation. In that 
case, the collision term is calculated as if the system were spatially homogeneous or as 
if the collisions could be treated as local. Then, the effect of spatial inhomogeneity is 
only retained in the advective (Vlasov) term which describes the evolution of the system 
due to mean- field effectfl This leads to the Vlasov-Landau-Poisson system which is 
the standard kinetic equation of stellar dynamics. To our knowledge, this equation 
has been first written and studied by Henon [23]. Henon also exploited the timescale 
separation between the dynamical time to and the relaxation time tn ^> to to derive a 
simplified kinetic equation for f(e,t), where e = v 2 /2 + $(r, t) is the individual energy, 
called the orbit-averaged-Fokker-Planck equation. In this approach, the distribution 
function /(r, v,t), averaged over a short timescale, is a steady state of the Vlasov 
equation of the form f(e, t) which slowly evolves in time (on a long timescale) due to the 
development of "collisions" (i.e. correlations due to finite N effects or graininess). Cohn 
[24] numerically solved this equation to describe the collisional evolution of star clusters 
and investigate the gravothermal catastrophe that was predicted by Antonov [25] and 
Lynden-Bell & Wood [26] on the basis of statistical mechanics. The local approximation, 
which is a crucial step in the kinetic theory, is supported by the stochastic approach of 
Chandrasekhar & von Neumann [10] showing the preponderance of the nearest neighbor. 
However, this remains a simplifying assumption which is not easily controllable. In 
particular, as we have already indicated, the local approximation leads to a logarithmic 
divergence at large scales that is difficult to remove. This divergence would not have 
occurred if full account of spatial inhomogeneity had been given since the start. On 
the other hand, the Fokker-Planck equation is based on a Markov assumption and this 

the Landau equation which conserves the energy (microcanonical description). 

|| This is not obvious on the form given by King, that is deduced from the work of Rosenbluth et al. 
[9]. However, the connection with the Landau equation can be made easily by using simple integrations 
by parts; see Ref. [18] for a recent discussion of that issue. 

% To our knowledge, the first explicit reference to the Landau equation appeared much later in the 
paper of Kandrup [19j . There is also a strange comment related to the work of Landau in the paper of 
Cohen et al. [TT] . 

+ The nowadays called "Vlasov equation" was introduced by Jeans [20] in stellar dynamics and by 
Vlasov |21j in plasma physics. See Henon [22| for interesting comments about the name that one 
should give to that equation. 
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assumption is not clearly justified for stellar systems. Memory effects can be important 
for self-gravitating systems because, as shown by Chandrasekhar [27], the temporal 
correlation function of the force decreases algebraically as 1/t instead of exponentially. 
This slow decay results in a logarithmic divergence of the diffusion coefficient for large 
times when one considers the Kubo formula. As shown by Lee [28], the spatial and 
temporal logarithmic divergences are equivalent, i.e. they are two manifestations of the 
same phenomenon. The effect of spatial and temporal derealization was investigated 
by Gilbert [29J, Severne & Haggerty [30], and Kandrup [19]. In particular, Kandrup 
derived a generalized Landau equation by using projection operator technics. Recently, 
Chavanis (3TJ [32] obtained this equation in a simpler manner from the BBGKY hierarchy 
or from a quasilinear theory as an expansion in 1/N in a proper thermodynamic limit. 
This generalized kinetic equation is interesting because it can take into account effects of 
spatial inhomogeneity and memory which are neglected in the previous approaches. It 
clearly shows which approximations are needed in order to recover the Landau equation. 
However, the generalized Landau equation remains extremely complicated for practical 
applications. 

Two-dimensional point vortices and non-neutral plasmas: It is well-known that 
a non-neutral plasma under a strong magnetic field is isomorphic to a system of 
point vortices in 2D hydrodynamics [33]. These systems are described by Hamiltonian 
equations in which the coordinates x and y of the particles are canonically conjugate [34] . 
The particles (charges or vortices) interact through a long-range potential that is played 
by the electric potential in a plasma or by the stream function in 2D hydrodynamics. 
The statistical mechanics of 2D point vortices was pioneered by Onsager [35J and further 
developed by Montgomery & Joyce [33] and Lundgren & Pointin [36] in a mean field 
approximation (recently, Eyink & Sreenivasan [37] have discovered that Onsager had 
first performed this mean field theory in unpublished notes). Kinetic theories of these 
systems were developed much later. Dubin & O'Neil [38] derived a kinetic equation 
for a 2D guiding center plasma by starting from the Klimontovich equation and using 
a quasilinear theory. Independently, Chavanis [39] derived a kinetic equation for 2D 
point vortices by using projection operator technics and, more recently [ID], by using a 
BBGKY-like hierarchy or a quasilinear theory. The kinetic equation derived by Dubin 
& O'Neil takes into account collective effects and can be considered as the counterpart 
of the Lenard-Balescu equation in plasma physics. The kinetic equation derived by 
Chavanis ignores collective effects and can be considered as the counterpart of the 
Landau equation in plasma physics and stellar dynamics. However, the approach of 
Chavanis is formulated for an arbitrary distribution of point vortices while the approach 
of Dubin & O'Neil is restricted to axisymmetric distributions. These kinetic equations 
are valid at the order 0(1/N) in a proper thermodynamic limit, so that they describe 
the evolution of the system on a timescale of order Ntn- For N — > +oo, we get the 
2D Euler equation which describes collisionless systems. At the order 1/N, the collision 
operator describes the development of correlations between point vortices. However, 
when we consider axisymmetric distributions, it is found that the collision operator 
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cancels out when the profile of angular velocity becomes monotonic. In that case, the 
evolution stops even if the system has not reached statistical equilibrium [391 El HQ] • 
This implies that the relaxation time of the system as a whole is larger than Ntp, where 
tn is the dynamical time, so that higher order correlations (three-body, four-body...) 
must be advocated [39]. It is also possible that the point vortex gas does not reach 
the Boltzmann distribution predicted by statistical mechanics. Indeed, the dynamics 
may be non-ergodic as discussed by Khanin [32]. On the other hand, using an analogy 
with stellar dynamics and Brownian theory, Chavanis [331 EHl 130] derived a Fokker- 
Planck equation describing the relaxation of a test vortex in a bath of field vortices at 
statistical equilibrium. The relaxation is due to a competition between a diffusion term 
and a drift term. The systematic drift [33], which is the counterpart of the dynamical 
friction [12] in stellar dynamics, is necessary to obtain the Boltzmann distribution at 
statistical equilibrium. The diffusion coefficient and the drift term are related to each 
other by an appropriate Einstein relation and they are inversely proportional to the local 
shear created by the field vortices. The properties of this Fokker-Planck equation have 
been studied by Chavanis & Lemou [31] ■ On the other hand, Chavanis & Sire [33] have 
studied the statistics of the velocity fluctuations arising from a random distribution of 
point vortices by using an approach similar to the one developed by Chandrasekhar & 
von Neumann [10] for the gravitational force, but leading to different results due to the 
lower dimension of space. The numerous analogies between the statistical mechanics 
and the kinetic theories of stellar systems and 2D vortices are discussed by Chavanis 
|45j*l . Other interesting kinetic theories of 2D point vortices have been developed by 
Nazarenko & Zakharov [36], Marmanis [31], Sire & Chavanis [3H] and Newton & Mezic 

HMF model: in the past, the kinetic theories of systems with long-range interactions 
have been essentially developed for plasmas, stellar systems and 2D vortices. Recently, 
there was a renewed interest from a part of the community of statistical mechanics for 
systems with long-range interactions. In that context, a simple model of systems with 
long-range interactions, the so-called Hamiltonian Mean Field (HMF) model, has been 
extensively studied [50]. It consists of iV particles moving on a ring and interacting 
via a cosine potential. It can be viewed as a one dimensional periodic system where 
the potential of interaction is truncated to one Fourier mode. This is certainly the 
simplest system with long-range interactions that we can imagine. However, it remains 
highly non trivial and, interestingly, it exhibits many features common to other systems 
with long-range interactions such as self-gravitating systems. The kinetic theory of the 
spatially homogeneous phase of the HMF model has been developed by Bouchet [51] . 
Bouchet & Dauxois [52] and Chavanis et at. [53] (Inagaki [53] first developed a kinetic 
theory of the HMF model but he arrived at incorrect conclusions). They considered an 

* As reported by Eyink & Srcenivasan [37] , Onsager first pointed out analogies between stellar systems 
and 2D vortices. In a Letter to Lin he wrote: "At negative temperatures, the appropriate statistical 
methods have analogues not in the theory of electrolytes, but in the statistics of stars...". However, 
Onsager did not develop this matter further. 
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expansion of the equations of the dynamics in powers of 1/N in a proper thermodynamic 
limit. For N — > +00, this leads to the Vlasov equation. At the order 0(1/N), this leads 
to the one dimensional version of the Lenard-Balescu equation. However, the Lenard- 
Balescu collision operator is known to vanish in one dimension [55| [56] . This implies 
that the relaxation time is larger than Nto, where to is the dynamical time. This 
result of kinetic theory is in agreement with direct numerical simulations that lead to 
£r ~ N 1 ' 7 tr, [57] or even t^ ~ e N t^ [58] depending on the initial conditions. These 
authors also considered the relaxation of a test particle in a bath of field particles. It 
is described by a Fokker-Planck equation, involving a diffusion and a friction, which is 
the one dimensional version of the Fokker-Planck equation derived by Hubbard [5] in 
plasma physics. It can also be deduced from the Lenard-Balescu equation by making 
a bath approximation^. Bouchet & Dauxois [52] used this Fokker-Planck equation to 
show that the temporal correlation function of the velocity decreases algebraically with 
time leading to (weak or strong) anomalous diffusion (related results were previously 
obtained by Marksteiner et al. [59] in the very different context of optical lattices). 
The theoretical results of Bouchet & Dauxois [52] are in very good agreement with 
direct numerical simulations in the situations considered [60]. Chavanis & Lemou [6"T] . 
developing the theory of Potapenko et al. [62] , used the Fokker-Planck equation to study 
the relaxation of the distribution function tail and show that it has a front structure 
moving very slowly (logarithmically) with time. Some features of the kinetic theory of 
the spatially inhomogeneous phase of the HMF model have been considered in Ref. [63J . 

A first objective of the present paper is to develop a unified kinetic theory of 
systems with long-range interactions. In particular, we show that, when collective 
effects are neglected, a generalized kinetic equation can be given for stellar systems 
[see Eq. (|2J)], for the HMF model [see Eq. (J2])] and for 2D point vortices [see Eq. ([52]i ]. 
This equation, which is not well-known, applies to spatially inhomogeneous systems (or 
non axisymmetric distributions of point vortices) and bears a clear physical meaning 
in terms of generalized Kubo relations. Furthermore, it avoids the use of Fourier- 
Laplace transforms and remains in physical space. For homogeneous distributions (or 
for axisymmetric distributions in the case of point vortices), it reduces to the Landau 
equation (or to its counterpart in 2D hydrodynamics). Therefore, our paper offers a 
new derivation of Landau-type equations from a simple formalism and develops the 
calculations in detail. A drawback of our approach is to neglect collective effects 
that would lead to Lenard-Balescu-type equations in the homogeneous case. A second 
objective of the paper is to use this general formalism to derive a kinetic equation, 
written in angle-action variables, describing spatially inhomogeneous systems with 

jj The Lenard-Balescu equation and the Fokker-Planck equation keep the same form in any dimension 
of space, and for any potential of interaction u(\r — r'Q, differing only in a multiplicative constant (2ir) d 
[18] . Therefore, they can be directly applied to the HMF model, and more generally to any system 
with long-range interactions. However, the explicit expressions of the diffusion coefficient and friction 
force depend on the dimension of space. Some useful expressions in d = 3, 2, 1 dimensions are given in 
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long-range interactions. This new derivation solves the shortcomings of our previous 
derivation [64] . We also derive general expressions for the diffusion and friction (or 
drift) coefficients of a test particle evolving in a bath of field particles and make 
contact with the expressions previously obtained in the literature. Finally, we use 
this kinetic theory to argue that the relaxation time is shorter for inhomogeneous (or 
high-dimensional) systems than for homogeneous (or low-dimensional) systems because 
there are potentially more resonances, and we compare this prediction with existing 
numerical results. For the HMF model, we argue that the relaxation time scales as N 
for inhomogeneous distributions and as e N for permanently homogeneous distributions. 
Phase space structures can reduce the relaxation time by creating some inhomogeneities 
and resonances. Similar results are expected for 2D point vortices. 



2. Stellar systems and HMF model 

2.1. The evolution of the system as a whole: the kinetic equation 

Let us consider an isolated system of particles with long-range interactions described by 
the Hamiltonian equations 

dvi dH dwi dH 

dt <9v,j ' dt dYi ' 

1 N 

/ i=l i<j 

where u(r — r') is a binary potential of interaction depending only on the absolute 
distance |r — r'| between the particles. We consider the proper thermodynamic limit 
N +00 in such a way that the normalized energy e = E/ (u*N 2 m 2 ) and the normalized 
temperature rj = (3Nm 2 u* are of order unity, where u* represents the typical value of 
the potential of interaction [65]. By a suitable normalization of the parameters, we 
can consider N — > +oo with V A ~£'/iV r ~/3~m~i and u* ~ 1/N. This is the 
usual Kac prescription which amounts to putting in front of the interaction energy 
(for rigorous mathematical results see (23 EI])- We can also consider N —> +oo with 
V^Er^fl/Nr^u*^! and m ~ 1/N. When collective effects are ignored, the 
kinetic equation describing the evolution of the distribution function f(r,v,t) at the 
order 0(1/N) is given by 

/( r v t _ T )£( r v t_ T ) ( 2 ) 
m 

where F(l — > 0) is the force (by unit of mass) created by particle 1 (with position r! 
and velocity vi) on particle (with position r and velocity v) and J-" M (1 — > 0) = F(l — > 
0) — (F)/N is the fluctuating force with respect to the average force (F)(r, t) = — V$(r, t) 
experienced by particle 0. On the other hand, G(t,t — r) is the Greenian constructed 
with the mean flow in phase space. The kinetic equation (T5]) is not well-known and this 
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is why we try to publicize it in this paper. It has been obtained by Kandrup [TH] from 
a projection operator formalism and by Chavanis [311 [32] from the BBGKY hierarchy 
and from a quasilinear theory. The BBGKY hierarchy is closed by neglecting the three 
body correlation function which is of order 0(1/N 2 ). The structure of this equation 
bears a clear physical meaning. It involves a diffusion term and a friction term. The 
coefficients of diffusion and friction are given by generalized Kubo formulae, i.e. the 
time integral of the temporal correlations of the fluctuating force. The collision term 
Cjv[/] (r.h.s.) is valid at the order 1/N so that it describes the "collisional" evolution of 
the system (ignoring collective effects) on a timescale of order Ntjj. For N — > +oo, we 
obtain the Vlasov equation in which collisions (more properly correlations) are neglected 
(for rigorous mathematical results see [68 | [69 | [70]). When coupled to an attractive long- 
range potential (e.g., the gravitational potential), the Vlasov equation can generate a 
process of violent collisionless relaxation on a few dynamical times to- A statistical 
theory of violent relaxation has been developed by Lynden-Bell [7_T] for stellar systems, 
but its domain of applicability is more general. Equation (j2|) is a non-Markovian integro- 
differential equation. It can describe the evolution of stellar systems, or the evolution 
of the HMF model, taking into account derealizations in space and time (i.e. spatial 
inhomogeneity and memory effects). If we make a Markovian approximation and extend 
the time integral to +oo, we obtain 



df_ df_ N-l,^ df d 
~dt +V '7h + ^J 

xG(t,t-r) 



J^(l^O)— + J^(0^1) ' 



^ .. .-^/(r.v.t^fr.v,,*). (3) 

The Markov approximation is not rigorously justified for self-gravitating systems because 
the force auto-correlation function decreases algebraically like 1/t [27], instead of 
exponentially. The Markov approximation is also expected to be incorrect for the 
HMF model and other systems with long-range interactions if we are close to the 
critical point. However, except for these situations, the Markovian approximation 
should be justified in the N — > +oo limit because the timescale ~ Nto on which 
f(r, v, t) changes is long compared to the timescale r corr for which the integrand in Eq. 
(|2]) has significant support. The Markovian equation <$5§ applies to possibly spatially 
inhomogeneous distribution functions, which is a specificity of systems with long-range 
interactions. Now, for spatially homogeneous distribution functions /(v, t), the kinetic 
equation reduces to 

df a 



X 



■er - --.^A^ut). (4) 

This is valid as long as the homogeneous distribution is Vlasov stable. In fact, the 
distribution f(v,t) will change due to the development of correlations and graininess 
effects (i.e. the r.h.s. of the kinetic equation) and, at some point of the evolution, the 
distribution may becomes Vlasov unstable. In that dynamical phase transition 
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from a homogenous state to an inhomogeneous state can take place. Such a transition 
has been illustrated numerically for the HMF model by Campa et al. |72J . Using the 
symmetry of the force F(0 — > 1) = — F(l — > 0), and the equations of the unperturbed 
trajectory 

v(t - r) = v(t) = v, r(t - t) = r(t) - v(t)r = r - vr, (5) 
corresponding to (F) = 0, the kinetic equation can be written 



(6) 



with 



r+oo r 

= J drj dnF"(l ^0,^(1 -»• 0,*-r). (7) 
Now, the force (by unit of mass) created by particle 1 on particle is given by 

F(l->0) = -m^(r-r 1 ). (8) 

OT 

Introducing the Fourier transform of the potential 

u(x) = / e^u(k) dk, u(k) = J e- ik " x u(x) ^L, (9) 

we get 

F(l ->■ 0, s) = -im J ke ik - (r(s) - ri(s)) w(k) dk. (10) 

Using the equations of motion (jSj), and introducing x = r — r x and w = v — Vi, we 
obtain 

F(l ^ 0,t - r) = -im |ke lk - (x ^ WT) n(k)dk. (11) 



Therefore 

r+oo 



= -™ 2 J q °° dr j dx j dk J dk ! k^ u e^^e-^'^u^uik!). (12) 
Using the identity 

w-f^W? < 13 > 

and integrating on x, then on k', we find that 

K i» = (2Ti) d m 2 J dr J rfkFre*- WT !i(k)ii(-k). (14) 

Since u(x) is real, we have u(—k) = u(k)*. Furthermore, since u{— x) = w(x), we have 
u(k)* = u(k). Therefore, u(k) is real and u(—k) = u(k). Therefore, 

f + OO f 

= (2vr)<Vi 2 J dr J dkPfcV kwT £(k) 2 . (15) 

Making the transformation r — > — r, then k — > — k, and adding the resulting expression 
to Eq. ( |T5|) . we get 

r+oo r 

= -{2Ti) d m 2 / dr / dk P£;V kwT w(k) 2 . (16) 
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Using the identity (1TB1) . we finally obtain 

K"" = -K(2n) d m 2 J dk k"k u 5(k ■ w)n(k) 2 . (17) 
Therefore, the kinetic equation takes the form 

% = 7r(27r) d m-^- / WW 50s. ■ w)n(k) 2 ( fi^J- - f^-) dv x dk, (18) 
dt v ; dW K 1 V ; V dv v J dv\) K 1 

where we have noted / = /(v, t) and f\ = /(vi,t). The calculation of K^ u is detailed 
in 



Appendix A[ Finally, the kinetic equation can be written 



df d f w 2 5^-wW f f df M\ 

m = Kd d^J ^ [f^ ~ f-Q^J d ^ (w) 

with 

r+ao r+oo 

K 3 = 8yr 5 m / k 3 it(k) 2 dk, K 2 = 8vr 3 m / k 2 u{k) 2 dk. (20) 
Jo Jo 

This is the original form given by Landau [2] for the Coulombian interaction (in d — 3). 

The Landau equation is derived here in a very different manner. This equation also 

applies to stellar systems if we make a local approximation. In that case, the distribution 

functions / = /(v, i) and f\ = /(vi,i) must be replaced by / = /(r, v,t) and 

/i = /(r, Vi,t) and we must restore the advective (Vlasov) term in the l.h.s. where the 

gravitational field is determined by the Poisson equation. Using (27r) 3 w(k) = —4irG/k 2 , 

this yields the Vlasov-Landau-Poisson equation 

df df ^ df ^ d fw 2 5P u - W Pw v f e df t df\\ , 
+ v ._L_ V $- 7 f = K— / (fJ-f.Jl)dv lt 



dt dr dv dvf 1 J w 3 V dv u dv", 

(21) 

A$ = 4ttG J /(r, v, t) dv, (22) 
where = 2-KmG 2 In A and In A = J +o ° * is the usual Coulomb logarithm that has to 



be regularized with appropriate cut-offs. It scales like In A ~ In A r (see, e.g., 

The kinetic equation generalizing the Landau Eq. fl20l by taking into account 
collective effects is the Lenard-Balescu equation [3j H] : 

dt 

where efk, w) is the dielectric function 



(23) 



r k • ^ 

e(k,o;) = 1 + (27r) d «(k) / -2s_ dv. (24) 

J cj — k ■ v 

The Landau equation is recovered by taking |e(k, k • v)| 2 = 1. In plasma physics, it has 
been shown that collective effects regularize the logarithmic divergence at large scales 
that appears in the Landau equation. The Landau and the Lenard-Balescu equations 
conserve mass M = f f drdv and energy (reducing to the kinetic energy E = f f^ drdv 
for a spatially homogeneous system) and monotonically increase the Boltzmann entropy 



Kinetic equations for systems with long-range interactions: a unified description 



12 



S = — f ±- In ^ drdv (.£/-theorem). The collisional evolution is due to a condition of 
resonance between the orbits of the particles. For spatially homogeneous systems, the 
condition of resonance encapsulated in the 5-function corresponds to k • v = k • vi with 
v 7^ Vi. For d > 1, the only stationary solution of the Landau and Lenard-Balescu 

2 

equations is the Maxwell distribution / = Ae~^ m ^ and these equations relax towards 
the Maxwell distribution for t — > +oo. Since the collision term is valid at the order 
0(1/N), the relaxation time scales like 

t R ~ Nt D , (d > 1). (25) 

This scaling, predicted in [73], has been confirmed for a two-dimensional plasma in 
[71| [T8] . For the Newtonian potential in d = 3 there is a well-known logarithmic 
divergence at large-scales in the Landau equation (|2"TT) . This divergence is due to the 
assumption of spatial homogeneity and to the linear trajectory approximation, and is 
not present in the generalized Landau equation (j2J). When this divergence is properly 
regularized, the relaxation time is found to be t R ~ (N/ hi N)tr> (see, e.g., [31]). This is 
the Chandrasekhar relaxation time. For one-dimensional systems, like the HMF model, 
the kinetic equation ( 1231) reduces to 

where we have used the identity 6(Xx) = r^6(x). Therefore, the collision term Cjv[/] 
vanishes at the order 1/N because there is no resonance. The kinetic equation reduces 
to df/dt = so that the distribution function does not evolve at all on a timescale 
~ Nto- This result is known for a long time in plasma physics and was discovered by 
Eldridge & Feix [55] (see also [SB]) for a one-dimensional plasma: "To first order in the 
small parameter g = {nD)' 1 , the detailed balance between drag and diffusion is valid not 
only at thermal equilibrium but for any stable function f(v). So the Maxwellization is, 
at least, a second-order effect in g and consequently a very slow process". This result has 
been rediscovered recently for the HMF model by Bouchet & Dauxois [52] and Chavanis 
et al. [53]. This implies that, for one-dimensional homogeneous systems, the relaxation 
time towards statistical equilibrium is larger than Nto- We therefore expect that 

t R >Nt D , (d=l). (27) 

For the HMF model, it has been found numerically that t R ~ N 1 ' 7 tu [S7J or t R ~ e N tr, 
[58] , depending on the type of initial conditions. For ID inhomogeneous distributions, 
the relaxation time can be reduced and approach the natural scaling Ntu associated to 
the generalized kinetic equation (see discussion in Sec. |5]). 

2.2. Test particle in a thermal bath: the Fokker-Planck equation 

We now consider a "test" particle (tagged particle) evolving in a steady distribution of 
"field" particles. The test particle has a stochastic motion. We assume the system to 
be spatially homogeneous. Let us call P(v, t) the probability density of finding the test 
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particle with velocity v at time t. The evolution of P(y,t) can be obtained from the 
Lenard-Balescu equation (123]) by considering that the distribution f\ of the field particles 
is fixed. Thus, we replace / = /(v, t) by P = P(v,t) and fi = /(vi,t) by fi = /(vi) 
where /(v) is any stable stationary solution of the Vlasov equation. This procedure 
transforms the integro-differential equation ( 1231) into the differential equation 

(25 



where e(k, w) is the dielectric function corresponding to the fixed distribution function 
/(v). Equation ( 128]) can be written in the form of a Fokker-Planck equation 

dP d ( OP \ 

D^T^-Prin, (29) 



dt dv^ \ dv u 
involving a diffusion term 

= 7r(2n) d m [ kWSQc ■ w) , . |2 A d Vl dk, (30) 

j |e(k, k-vjp 

and a friction term due to the polarization 

if = K* = <^) d ™ J *rV6(k • w) |e( ff )2 v)|2 || dv^fe. (31) 

The diffusion term can be directly derived from the Kubo formula, and the friction 
force due to the polarization can be directly obtained from a linear response theory, as 
shown by Kandrup "~~5] and Chavanis [21] • Since the diffusion coefficient depends on 
the velocity v of the test particle, it is useful to rewrite Eq. ( 129]) in a form that is fully 
consistent with the general Fokker-Planck equation 

dP l d 2 ( p {Av»Av») \ d f p {^)\ (32) 



dt 2dv^dv u \ At J dv»\ At 
with 

{Av^Av v ) _ (A^> _ dD>» 

2At ~ ' At ~ 9u» +V friction- W 

Substituting Eqs. (130]) and (13*1]) into Eq. (|33j) . and using an integration by parts, we 
find that the diffusion and friction coefficients are given by 

<A ^ = ^ / k"kS(k ■ w) |e(k fi( k k)2 v)|2 /i dv.dk, (34) 

The two forms (129]) and (132]) of the Fokker-Planck equation have their own interest. 
The expression (132]) where the diffusion coefficient is placed after the two derivatives 
d 2 (DP) involves the total friction force Ff r i Ction = (Av 11 ) / At and the expression (129]) 
where the diffusion coefficient is placed between the derivatives dDdP isolates the part 
of the friction rf = F^ ol due to the polarization [32] . 
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In stellar dynamics, the coefficients of diffusion and friction were first obtained by 
Chandrasekhar 0, [7J by considering the mean effect of a succession of binary encounters. 
They were also computed by Cohen et al. [H] and Rosenbluth et al. [9] in a plasma 
physics context by neglecting collective effects. This leads to expressions that are 
consistent with Eqs. ( )34l) and (1351) with |e(k, k ■ v)| 2 = 1 (see, e.g., [18] for a discussion 
of this link). These expressions can also be obtained directly from the Hamiltonian 
equations of motion by making a systematic expansion of the trajectory of the particles 
in powers of 1/N in the limit N — > +oo as shown in Appendix A of [32]. On the other 
hand, the expressions (1341) and (I35I) for the diffusion and the friction, properly taking 
into account collective effects, were first derived by Hubbard [5] in plasma physics by 
directly evaluating (Av^Av u ) and (Aw M ). Collective effects encapsulated in the dielectric 
function remove the logarithmic divergence that otherwise appears at large scales. The 
work of Hubbard [5] was done independently from the works of Lenard [3] and Balescu 
[I] that were published at about the same time, but they are clearly connected since 
the equations ( IMl) and ( 1351) of Hubbard [5] can be derived from the Lenard-Balescu 
equation, and vis et versa, as explained above (see also discussion in [32]). 

Let us now consider particular cases. 

• In d = 1, the expressions of the diffusion and friction simplify into 



ku{k) 



2 



«=^»1 WMr dk - (36) 

<M = C>) + DW ^, (38) 

and the Fokker-Planck equation fl28|) reduces to 
OP _ d 
dt dv 



\dv dv ) \ 

The diffusion and the friction coefficients were derived by Eldridge Sz Feix [55] for one 
dimensional plasmas and by Bouchet & Dauxois [52] for the HMF model by directly 
evaluating {(Av) 2 ) and (Av) (see [50] for detailed calculations). Alternatively, Chavanis 
[73] (see also [531 EH]) obtains Eqs. ( |36|) -( 139|) as a particular case of the general 
expressions (I28l -(l35"l) valid in any dimension of space. 

• If we neglect collective effects, the diffusion and friction coefficients reduce to 

D» v = Ti{2ix) d m J k^k u 5(k ■ w)w(k) 2 /i d^dk, (40) 
jf = F; ol = ix{2n) d m J k^k u 5{k ■ w )^(k) 2 |^ d^dk, (41) 

Faction = <1n) d m J kWfr (A _ JL\ 5 {k ■ w)u(k) 2 dv x dk. (42) 
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After a series of elementary transformations, we obtain 
8D^ V r 8 

— = Ti(2n) d m J k^—6(k • w)£(k) 2 /i dv x dk, 

= -n(27r) d m [ k"k u — 5(k ■ w)u(k) 2 /i dvidk, 

J OVi 

= 7r(27r) d m J Wk u 5(k ■ w)u(k) 2 ^ d^dk = 77". (43) 
Combining Eq. ( 1331) with Eq. (1431) . we find that 

^ = V, i-e. F At-ta = 2F^ (44) 

so that the friction force F friction is equal to twice the friction due to the polarization 
F po /. This explains the difference of factor 2 in the calculations of Chandrasekhar [7] 
and Kandrup [19] (see the anomaly noted by Kandrup [IH] p. 446). 

• For a thermal bath, the field particles have the Maxwell distribution of statistical 
equilibrium 



■2 



/(Vl) = (¥) d/2pe "^- (45) 



Substituting 

dfi 



-Pmftf, (46) 



8v\ 

in Eq. ( 13~TT) . we obtain 

^ = -Pn(2Tx) d m 2 J F(k • Vl )5 [k ■ (v - Vl )] | 6(k M ^ )2 v)[2 /i ^idk. (47) 
Using the ^-function to replace k • vi by k • v, we get 

if = -f3ii(2n) d m 2 v v [ k»k u 5 [k ■ (v - Vi)] , tt^t^ rpr /i dv^k. (48) 
J |e(k,k-v)p 

Comparing the resulting expression with Eq. ( 1301) . we find that 

^ = -fimD^v", (49) 

which is the appropriate Einstein relation for our problem. Note that it is valid for the 
friction force due to the polarization F po i, not for the total friction force Ff riction . The 
Fokker-Planck equation (128]) takes the form 



8P 8 



(50) 



8t dv 11 

where the diffusion coefficient is given by Eq. ( 1301) with Eq. ( 1451) . This is similar to the 



Kramers equation except that the diffusion coefficient is a tensor and that it depends 
on the velocity. Some simplified forms of the diffusion coefficient are given in [18] in 
different dimensions of space d. 
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3. Two-dimensional point vortices 

3.1. The evolution of the system as a whole: the kinetic equation 

We consider a system of point vortices described by the Hamiltonian equations 
dxi dH dyi dH 

^ dt dyi ' ^ dt d%i ' 

= 7 2 EMr l -r J ). (51) 

i<j 

The usual interaction corresponds to «(r - r') = — =- In |r — r'| but we let u arbitrary 
in order to treat more general cases. We consider the proper thermodynamic limit 
N — > +oo in such a way that V~i?~/3/A r ~l and 7 ~ 1/N [65] (for rigorous 
mathematical results see [761 EH US CHI EH])- When collective effects are ignored, the 
kinetic equation describing the evolution of the smooth vorticity profile u(r,t) at the 
order 0(1/N) can be written 



x 



^(i^o)A + v ^i)A 



CO 



u(r,t-r)-(r u t-r), (52) 

7 



where V(l — > 0) is the velocity created by point vortex 1 (located in ri) on point 
vortex (located in r) and V M (1 — > 0) = V(l — >■ 0) — (V)/N is the fluctuating velocity 
with respect to the average velocity (V)(r, t) = — z x V^(r,t) at the location of point 
vortex 0. On the other hand, G(t,t — r) is the Greenian constructed with the mean 
flow. This kinetic equation has been obtained by Chavanis [391 SO] from a projection 
operator formalism, a BBGKY-like hierarchy and a quasilinear theory. The BBGKY 
hierarchy is closed by neglecting the three body correlation function which is of order 
0(1/N 2 ). The structure of this equation bears a clear physical meaning. It involves a 
diffusion term and a drift term. Furthermore, the coefficients of diffusion and drift are 
given by generalized Kubo formulae, i.e. the time integral of the temporal correlations 
of the fluctuating velocity. The collision term Cat[u/| (r.h.s.) is valid at the order 1/N so 
that it describes the "collisional" evolution of the point vortex gas (ignoring collective 
effects) on a timescale of order Ntr>. For N —> +00, we obtain the 2D Euler equation 
in which collisions (more properly correlations) are neglected. The 2D Euler-Poisson 
system can generate a process of violent relaxation towards a quasi stationary state 
as described by Miller [81] and Robert & Sommeria [82]. This is the hydrodynamical 
counterpart of the Lynden-Bell [7_1] theory of violent relaxation in stellar dynamics (see, 
e.g. [831 HH], for a description of this analogy). Equation ( 152]) is a non Markovian 
integro-differential equation. It can describe the dynamics of point vortices taking into 
account delocalizations in space and time (i.e. non-axisymmetry and memory effects). 
This kinetic equation is the vortex analogue of the generalized Landau equation ([2J) for 
stellar systems. If we make the Markovian approximation and extend the time integral 
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to +00, we obtain 



dcu N — 
~dt + ~W 



i(V> • £ 

or 



d 



dr / dnV^l ->0)G(M-r) 



^(i^o)A + v^(o^i)A 



u(r,t)-(r 1 ,t). 

7 



(53) 



The Markovian approximation may not be justified in any situation since it has 
been shown numerically that point vortices can exhibit long jumps (Levy flights) and 
correlations (see Sec. [5] for a more detailed discussion). However, for N — > +00, the 
Markovian approximation should be justified because the timescale ~ Nto on which 
cu(r, t) changes is long compared to the timescale r corr for which the integrand in Eq. 
( 152"j) has significant support. The Markovian equation ( 153"j) still applies to possibly non- 
axisymmetric distributions, which is the norm for 2D flows. Now, for axisymmetric 
distributions u(r,t), the kinetic equation reduces to 
du 1 d r+00 
~dt 



r dr j 



r+OO r 

J dr J dnV r (l 0)G(t, t - r) 



K(i^o)— + v ri (o^i)— 

or ori 



u(r, t)-(r h t). 
7 



(54) 



Using V ri (0 — > 1) = —^V r (l —> 0) (see below) and the equation of the unperturbed 
trajectory 

r(t-r) =r(t) = r, 9(t - t) = 9(t) - Q(r, t)r = 9 - fi(r,t)r, (55) 

corresponding to (V) = V(r,t)eg with V(r,t) = Q(r,t)r (where Q is the angular 
velocity), the kinetic equation can be written 



dU - 1 9 [ +0 ° r dr K t 1 d 
dt r dr Jo 1 1 \r dr 



1 d \ . s cu . . 
r ± or 1 J 7 



with 



K = r 2 dr d9 1 V r (1^0,t)V r (1^0,t-r) 
Jo Jo 



(56) 



(57) 



In consistency with the Markov approximation, we have assumed that the angular 
velocity f2(r, t) does not change on the timescale r corr on which the velocity correlations 
have their support. The velocity created by point vortex 1 on point vortex is given by 
V(l — >• 0) = — 7Z x Vw(r — ri) leading to 

7 du 



K(l->0) = - — (r 
r 09 



(5? 



The potential of interaction is of the form u(r — ri) = u(r,ri,9 — 61). It depends on 
cos(6* — 61) and is invariant with respect to the interchange of and 1. Introducing the 
Fourier transform with respect to the angles 



u(r, n, 



2n 

u m [r,r 1 ) = —i u(r,ri,(j))cos(m(, 
2ir Jo 



(59) 
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we get 

V r {\ -)> 0,s) = i^Y, memieis) ~ 9lis)) ^(r(s),r 1 (s)). (60) 



rls) 



Using the equations of motion ( )55l) . and introducing <fi = 9 — 6\ and AQ = Q(r,t) — 
Q(ri, t), we obtain 

V r (l 0, t - r) = i 1 - me im ^ An ^u m (r, n). (61) 

r m 

Therefore 

r+oo /■2tt 

K = - 1 2 dr #Vmrie im< ^e in ^ Anr ^ m (r,ri)w n (r,ri). (62) 
Using the identity 

Sn,o= / e^f , (63) 

JO Z7T 

and integrating on 0, then summing on n, we find that 

K = 27T 7 2 / dTY.™ 2 e imAnT u m (r 1 r 1 )u_ m (r 1 r l ). (64) 

Since u_ m (r, 77) = M m (r, ri), the foregoing equation can be rewritten 

r+oo 

K = 2tt 7 2 / dr ]T m 2 e imAnr u m {r, nf. (65) 
Jo m 

Making the transformation r — > — t, then m — > —m, and adding the resulting expression 
to Eq. (|65|) . we get 

/+oo 
c/r^mV^^^n) 2 . (66) 
-oo m 

Using the identity ([TBI , we finally obtain 

K = 2tt 2 7 2 m 2 5{mAQ)u m {r, n) 2 . (67) 

m 

Therefore, the kinetic equation takes the form 

du 2 1 9 r+oo A ^ i 

_ = 2vr 7-— / n cfrix(r, ri)d(i2 - ih) -wi- w-— 

ar r or Jo \r or n or\ y 



(6? 



with 

+oo 

x( r , r i) = M^™( r ; r i) 2 , (69) 



m=—oo 



where we have used the identity S(Xx) = ^j5(x) and noted oj = u(r,t), U\ = u(r^t), 
Q = Q(r,t) and Qi = Q(ri,t). The calculation of xi r i r i) is detailed in Appendix B for 
the ordinary potential u(r — r') = — ^-ln(r — r'). Equation (j68p has been derived by 
Chavanis [39]. It is the vortex analogue of the Landau equation fTTBl . A more general 
kinetic equation, taking into account collective effects, has been derived by Dubin & 
O'Neil [HE] by a different method. It is the vortex analogue of the Lenard-Balescu 
equation (123]) . The kinetic equation ( 168]) conserves circulation T = fudr, energy 
E — | J uip dr and angular momentum L = J ur 2 dr, and monotonically increases the 
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Boltzmann entropy S — — f - In - dr (if -theorem). The collisional evolution is due to 
a condition of resonance between distant orbits of the point vortices. For axisymmetric 
systems, the condition of resonance encapsulated in the ^-function corresponds to 
Q(r, t) = Q(ri,t) with r ^ r\. The evolution stops when the profile of angular velocity 
becomes monotonic (so that there is no resonance) even if the system has not reached 
the statistical equilibrium state given by the Boltzmann distribution [33]. This "kinetic 
blocking" has been illustrated numerically in [JT]. Indeed, the Boltzmann distribution 
uj = Ae~ l3 ' y ^ + ^ nLr2 ' > is not the unique steady state of the kinetic equation ( 1681) : any 
vorticity distribution associated with a monotonic profile of angular velocity is a steady 
state of Eq. (1681) . In particular, if the profile of angular velocity is initially monotoni df f . 
the collision term CVfw] vanishes at the order 1/N because there is no resonance [39J. 
The kinetic equation reduces to dcu/dt = so that the vorticity does not evolve at all on 
a timescale ~ Ntn- This implies that, for axisymmetric distributions of point vortices, 
the relaxation time towards statistical equilibrium is larger than Nt^- We therefore 
expect that 

tn > Nto (axisymmetric flows). (70) 

For non axisymmetric distributions, the relaxation time can be reduced and approach 
the natural scaling Nto associated with the generalized kinetic equation ( 1521) (see 
discussion in Sec. [5]). In any case, for N 3> 1, the relaxation towards the Boltzmann 
distribution is a very slow process. The possible slow timescale of mixing was pointed 
out by Onsager in a letter to Lin [37]: "I still have to find out whether the processes 
anticipated by these considerations are rapid enough to play a dominant role in the 
evolution. " and it is now confirmed by the kinetic theory (as far as we know, the 
process of violent relaxation [7T| |8T| 182] was not foreseen by Onsager who focused on 
the ordinary statistical equilibrium state). In fact, the scaling of the relaxation time with 
N is still not known for axisymmetric distributions and it would be interesting to obtain 
it numerically. It is not even clear whether the system of point vortices truly relaxes 
towards the Boltzmann distribution predicted by statistical mechanics |33j . Indeed, as 
discussed by Khanin [12], the dynamics may be non-ergodic. This was also a concern 
of Onsager: "We inquire about the ergodic motion of the system, " Onsager wrote to Lin 
[37]. 

3.2. Relaxation of a test vortex in a thermal bath: the Fokker-Planck equation 

We now consider the relaxation of a "test" vortex (tagged particle) evolving in a steady 
distribution of "field" vortices. The test vortex has a stochastic motion. Let us call 
P(r,t) the probability density of finding the test vortex at position r at time t. For 
axisymmetric distributions, the evolution of P(r,t) can be obtained from the kinetic 
equation ( 1681) by considering that the distribution of the field vortices u\ is fixed. Thus, 
we replace u = ou(r, t) by P = P(r, t) and U\ = w(ri, t) by u± = u(ri) where u(r) is any 

ff This is the generic case because the process of violent relaxation usually leads to an axisymmetric 
QSS with a monotonic profile of angular velocity 
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stable stationary solution of the 2D Euler equation. This procedure 
the integro-differential equation f|68|) into the differential equation 



dP 
~dt 



2?r 2 7 



ld_ 

r dr 



+00 



ridrix{r, r\)5(Q 




ri dri 



transforms 

P(r,t)u(n). 
(71) 



This is the vortex analogue of Eq 
Fokker-Planck equation 

dP _ 1 d 

dt r dr 
involving a diffusion term 

2tt 2 7 



281) . Equation ( I7TT) can be written in the form of a 



dP 

r[ D— 

or 



Pr] 



D 



r!drix{r,ri)5(Q - Qi)cj(ri), 



and a drift term due to the polarization 

2tt 2 7 



V 



y P oi = 2^7 f + ™ driX ^ ri)6 ( Q _ Qi) dU( ri) 
r Jo dr 



(72) 



(73) 



(74) 



The diffusion coefficient can be directly derived from the Kubo formula and the drift 
velocity due to the polarization can be directly derived from a linear response theory 
[4~3l 1391 SO] . Since the diffusion coefficient depends on the position r of the test vortex, 
it is useful to rewrite Eq. flTTj) in a form that is fully consistent with the general Fokker- 
Planck equation 



dP 
~dt 



]_d_ 

2r dr 



dr \ At 



1 d 



r dr\ 



frP<£> 



At 



with 



<(Ar) 2 > 



(Ar) dD 



dr ^ ^ 



ydrift 



(75) 



(76) 



2At ' At 

Substituting Eqs. ( 1T3]) and (1741 in Eq. f!76|) . and using an integration by parts, we find 
that the diffusion and drift coefficients are given by 

((Ar) 2 ) _ 27^7 H 
r 2 Jo 



2At 

((Ar)) 



ridrix(r,ri)5(fi - fii)w(ri), 



(77) 



2vr 2 7 / 
Jo 



fid Id 

rrxdriui [-- 

\ r dr 



1 



X^r^n-n,)-. (78) 



At Jo yr9r r\dr\j 1- 

The two forms (17"2]) and (175|) of the Fokker-Planck equation are interesting to consider 
as they involve respectively the drift due to the polarization V^° l or the total drift V^ 1 " 1 * 1 
[4"0] . Expressions (I77|) and flTHj) for the diffusion coefficient and the drift term can be 
obtained directly from the Hamiltonian equations, by making a systematic expansion of 
the trajectories of the point vortices in powers of 1/N in the limit N — > +00 as shown 
in Appendix C of [4"T] . 

Let us now consider particular cases. 
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• If the profile of angular velocity of the field vortices Q(r) is monotonic, using the 
identity 5(Q — Q\) = 5(r — ri)/\Q'(r)\, we find that the expressions of the diffusion and 
drift simplify into 

D(r)=27r 2 7 |^(r), (79) 

yfr, r) ,, „ „. .din a; . . 

^ = ^ 2 7f^y|^) = ^(r)^-, (80) 

(Ar) . . . „. . din a; , . 

i_z = zy(r) + D(r) _ (81) 

where E(r) = rQ'(r) is the local shear. The Fokker-Planck equation (ITTj) can then be 
written 



rD(r)[ — P— — 

V or dr 



(82) 



dP _ 1 d 

dt r dr 

with a diffusion coefficient given by Eq. (1791) . Equation ( 182|) is the vortex analogue of 
Eq. ( |39l) . The expressions ( 179]) and (180]) of the diffusion coefficient and of the drift, t/iat 
are inversely proportional to the shear, were first derived for a thermal bath in [33] and 
for an arbitrary distribution of the field vortices in [39] (see remark before Eq. (123)). 
For the usual potential of interaction (see Appendix B ), we have x( r ? r ) = §~2 h 1 A where 



In A = Z)m=i * s a Coulombian logarithm that has to be regularized appropriately. It 
scales like In A ~ |hi A" (see, e.g. [4"T]). 

• For a thermal bath, the field vortices have the Boltzmann distribution of statistical 
equilibrium 

u{n) = A ie -^' {ri \ (83) 

where ip'(ri) = ip(ri) + ^rf is the relative stream function taking into account the 
invariance by rotation [41]. We have 



^ = -^(n)^ = ^WCnCn) - n^n, (84) 

where we have used Q(r) = V(r)/r = —(l/r)dip/dr. Substituting this relation in Eq. 
( 1741) . we obtain 

r Jo 

Using the 5-function to replace Q\ by f2, then using Q(r) — = —(l/r)dip'/dr, and 
comparing the resulting expression with Eq. (!73|) . we finally find that 

T) = -Dh^- (86) 

The drift is perpendicular to the relative mean field velocity (V) = —(dip' /dr)eg and 
the drift coefficient satisfies an Einstein relation £ = D/3^y [12] • We stress that the 
Einstein relation is valid for the drift V^° l = r) due to the polarization [3D], not for the 
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the Fokker-Planck equation (1721) 



dib 1 

\ or dr 



7) 



total drift V^ ri ^. For a thermal bath, using Eq. 
takes the form 

dP 1 d r ,/0P 

= fj_)ir\\ 

dt r dr [ 

where D(r) is given by Eq. (1731) with Eq. ( |83|) . Of course, if the profile of angular 
velocity of the Boltzmann distribution is monotonic, we find that Eq. (180j) with Eq. 
(183|) returns Eq. fl86l) with a diffusion coefficient given by Eq. f!79|) with Eq. fl83|) . Note 
that the systematic drift V po i = —Dft'y'Vip' of a point vortex [13] is the counterpart of 
the dynamical friction F po i = —D\\j3m\ of a star [7J and the Fokker-Planck equation 
(J87l) is the counterpart of the Fokker-Planck equation (J50l) . 

We have presented the results for axisymmetric distributions, but similar results 
can be obtained for unidirectional flows H31 [391. 



4. Kinetic equation with angle-action variables 

4-1. Evolution of the system as a whole: the kinetic equation 

The kinetic equation (TT8|) derived in Sec. [2] is valid for spatially homogeneous systems. 
A manner to describe spatially inhomogeneous systems is to use angle-action variables 
(0, J). A kinetic equation for f(J,t) has been derived by Chavanis [61] by using two 
different strategies: (i) by incorporating in a Fokker-Planck equation the expressions of 
the coefficients of diffusion and friction obtained by Valageas [El] or, (ii) by starting from 
the Klimontovich equation, using a quasilinear theory, and solving the equations with 
Laplace-Fourier transforms. However, the second derivation presents some shortcomings 
that were explicitly discussed in [61]. In particular, a factorization hypothesis is used 
without precise justification. Here, we show that, when collective effects are neglected, 
the kinetic equation in angle-action variables can be derived from the general equation 
df d f+°° 



dt dJ 

x 



r+oo f 

J dr J # lC ZJiF(l Q)G(t, t - t) 



F(l4 0)^- + f(0 4l) ' 



dJ y 'dJ x 



/(J,t)^(Ji,t), 
m 



with 



3n 

F(l -+ 0) = -m-Qr(x(J, 0) - x(J u 00). (89) 

Not only the present approach solves the problems of our former approach |64j . but it 
also shows a nice unity with the kinetic theories developed in Sees. [2] and [3] Indeed, 
the kinetic equation (1581 is the analogue of the kinetic equations (El) and ( 1531) . For 
simplicity, we consider one- dimensional systems but we shall indicate at the end how 
our results can be generalized in d dimensions. Introducing the Fourier transform of the 
potential with respect to the angles 

u{x{J,<f>)-x{J u <f> 1 ))= Y. A nM JOe^-^, (90) 

71, ni 
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where 

(91) 



we get 



An, ni (J,Ji) = 7333/ «(x(J,0) -x(J 1 ,0 1 ))e-^-^ 1 ) 

n,ni 

Substituting this expression in the kinetic equation (I88p . we obtain 
df 9 9 



i- 



(92) 



dt dJ Jo J n ni Mi 



A lA (J, Ji)lS l +- h +tjj + A lul (J u J)he 



i(h<t>i-l<t>) 



d_ 
d~h 



/(J,t)^(J 1; t). 
m 



(93) 



With angle-action variables, the equations of motion take a very simple form 

J(t -t) = J{t) = J, c(>(t -t) = <j){t) - fi(J, t)r = <j> - fi(J, t)r, (94) 

where 27rf2 _1 (J, £) is the period of the orbit with action J. Substituting these relations 
in Eq. (193]) and making the transformations / — > —I and li — > —l\ in the second integral 
(friction term), we obtain successively 

r) /"+00 r 

/ dr / dfadh E 4vu Ji)ne^-^ 

n.rai 



5/ 

— - = — m - 



xe 



i(Ut>(t-T)-h(f>l(t-T)) 



d d 
AhU Jt)l^j - A_ h ^(J u J)h-Q^ 



f(J,t)t(J h t), 
m 

(95) 



and 



at dJ Jo J nni ^ 



-i(m(j,t)-Zin(Ji,t))7 



d d 
A h (J, J x )l-Qj - A- h -i{J x , J)h-Qj-^ 



It is easy to establish that 

(Ji, J) = A 

—n,—n\ 

(J, Ji) — A ntTll (J, Jl)*. 
Therefore, the kinetic equation can be rewritten 



f(J,t)J-{J u t). 
m 

(96) 
(97) 



df 

m=- m dl 



d r+°° 



/*+oo r 

/ dr / Jx 2 £ A n , ni (J, J^n^^+O^-^^+'i)^ 



9 , 9 



/ 



Xe - W ,*)-^(^))r AMi(J) J i} Z - _ h - /(J >t )i_(J ljt ). (98) 



dJ d.h 



m 
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Integrating on (pi and using the fact that average distribution does not depend on 0, we 
get 



df 9 2 d 
-— = 2irm 
dt 



Q r+oo r 

— dr / dJi V \A nni (J, J 1 )\ 2 n 



Xe i(nn(j, t )- ni a(J 1 ,*))r / - _ ni ±\ f ^ t) L {J ^ t y (99) 




d J i J m 

Making the transformation r — > — r, then (n, n\) — > (— n, — n\), and adding the resulting 
expression to Eq. (199 p . we obtain 

~c% = 711112 t^J f dT f dJl ^ l^ n ' ni ^' Jl )\ 2n 




Xe i(nn(J,t)-n 1 n(J ll t))r ( _ ^ _^_\ /( J } J 1} f). (100) 

Using the identity (fl3l) . we finally obtain the kinetic equation 



5/ 9 2 d 



f dJ x J2 \ A n, ni (J, Ji)\ 2 n 

n,ni 



x5(nQ(J,t) - mfi(Ji,t)) - , (101) 



which returns the result obtained in [61]. However, it is derived here in a more 
satisfactory manner. In d dimensions, n, J, Q are replaced by vectors and the coefficient 
2tt 2 is replaced by 7i(2n) d . 

The kinetic equation f II 1 j) conserves mass and energy and monotonically increases 
the Boltzmann entropy (if-theorem) [M]- The collisional evolution is due to a condition 
of resonance encapsulated in the 5-function. The evolution stops when this condition 
of resonance cannot be satisfied anymore even if the system has not reached the 
statistical equilibrium state given by the Boltzmann distribution. Indeed, the Boltzmann 
distribution is not the only steady state of Eq. f II 1 j) : any distribution satisfying 
nQ(J,t) 7^ nifl(Ji,t) for any couples (n, J) ^ (m, Ji) is a steady state. However, 
the important point that we want to make here (it will be discussed further in Sec. |5]) is 
that, for inhomogeneous systems, there generically exists much more resonances than for 
homogeneous systems. In particular, for one dimensional homogeneous systems, there 
is no resonance at all (see Sec. [2]). We therefore expect that the relaxation time will 
be reduced for inhomogeneous ID systems as compared to homogeneous ID systems. 
In fact, we cannot conclude that Eq. (llOip tends to the Boltzmann distribution (since 
it is not the only steady state) but since entropy increases while energy and mass are 
conserved the system is expected to approach the Boltzmann distribution on the natural 
timescale Ntu on which Eq. (jlOip is valid, provided that there are enough resonances. 
Indeed, due to "collisions" and resonances, the system becomes "more mixed". These 
arguments will be further discussed in Sec. [5] in the light of existing numerical results. 
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4-2. Relaxation of a particle in a thermal bath: the Fokker- Planck equation 

Implementing a test particle approach like in Sees. !2.2l and [3T2l we find that the equation 
for P(J,t) is 

OP d r 

— = 2vr 2 m— J dJ x £ J i)\ 2n 

x6(nn(J) - mfi(Ji)) Ufi^r - n i p jJ^j • (102) 

Equation (11021) can be written in the form of a Fokker-Planck equation 

dP d ( dP \ . , 

W = aj( D aJ- p »y (103 » 

involving a diffusion term 

D = 27t 2 m J dJ x £ K, ni (J, J^nH^J) - mn(Ji))/i, (104) 



n,ni 



and a friction term due to polarization 

V = F po1 = 27r 2 m f dJ x £ \A^ ni (J,Ji)\ 2 nS(nQ(J) - n^J^n^. 

J n,ni dJi 

(105) 

The ordinary Fokker-Planck equation is 

dp i # f ((Ajy) \ s / p (aj}\ 

sT = 2 a/H^^ p J " a7 l p ^rJ ' (106) 

with 

Substituting Eqs. (I104p and (I105[) in Eq. (11071) and using an integration by parts, we 
find that the diffusion and friction coefficients are given by 

((AJ) 2 ) 



27r 2 m / dJxh J2 \A n , ni (J, Ji)\ 2 n 2 5(nQ(J) - rnfi(Ji)), (108) 

J n,ni 

^ = 2vr 2 m J dJxh £ - m-^-) \A n , ni (J, J^nd^J) - mfi(Ji)). 

(109) 



These expressions can be obtained directly from the Hamiltonian equations of motion 
by making a systematic expansion of the trajectory of the particles in powers of 1/N in 
the limit N -> +oo [83]. 

Let us consider particular cases. 

• If we consider a distribution of field particles /( J\) such that there is no resonance: 
nfl(J) 7^ nifl(Ji) for any couple (n, J) ^ (n l5 J x ), we first get 

D = 2-K 2 m [ dJ 1 ^|A nin (J,J 1 )|V5(n(fi(J)-fi(J 1 )))/ 1 . (110) 
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Then, using the identity 5(n(Q — f^)) = 8(J — J\)/\nQ!{J)\, we find that 

\^-n,n{Jj J)\ 2 



D(J) = 2ir 2 mf(J)J2 



n 



Similarly 



V = 2ir 2 mf'(J) £ 



n\ 



|0'(J)| 

A nt n(J, J)\ 2 



Using Eq. (I107p . we obtain 
(AJ) 



D'(J) + D(J) 



|fi'(J)| 
dlnf 



D{J) 



dhxf 
dJ ' 



At v ' dJ 

Finally, the Fokker-Planck equation f 1 1 U 3 j) can be written 

dP _ d 
~dt~dl 



D(J) 



OP dlnf 
P - 

dJ dJ 



[III] 



(112) 



(113) 



(114) 



with a diffusion coefficient given by Eq. (IllOp . 

• If the field particles are at statistical equilibrium (thermal bath), their distribution 
is the Boltzmann distribution 



/(Ji) = Ae-P me{Jl \ 



(115) 



where e(J) is the individual energy of the orbit with action J. Using de/dJ = fi(J), we 
find that 



4tL = -Pmf 1 %L = -/3mfiCl(J 1 ). 
dJi dJ\ 



(116) 



Substituting this relation in Eq. f)105p . we obtain 

n = -27r 2 m 2 f3 f dJ x £ \A n , ni (J, Jx)| 2 n5(nfi(J) - n 1 ^(J 1 ))n 1 / 1 ^(J 1 ). 

(117) 

Using the 5-function to replace nifl( Ji) by nQ( J) and comparing the resulting expression 
with Eq. fTTTHj) . we finally get 

r] = —Df3mQ(J), (118) 

which is the appropriate Einstein relation for our problem. For a thermal bath, using 
Eq. flllSp . the Fokker-Planck equation (11031) can be written 

dP _ d 
~dt ~dJ 

where D(J) is given by Eq. flTTMl) with Eq. ffTT5l) . 



BP \ 
D(J)( — + PmP{l(J)) 



(119) 



5. Discussion and conclusion 



The kinetic theory of systems with long-range interactions leads to the following 
scenario: 

1. The first stage of the dynamical evolution is described by the Vlasov equation. 
Due to violent relaxation and phase mixing, the coarse-grained distribution function 
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/(r, v, t) converges towards a steady state /(r, v) which is a dynamically stable (robust) 
steady state of the Vlasov equation. The convergence takes a few dynamical times 
to independent of N. The fine-grained distribution /(r, v, £) develops filaments at 
smaller and smaller scales and does not achieve a steady state (presumably). This 
collisionless relaxation explains the rapid emergence of coherent structures in stellar 
dynamics (galaxies), two-dimensional turbulence (jets and vortices) and in the HMF 
model (quasistationary states). 

2. If the system mixes efficiently, the QSS is given by the Lynden-Bell statistical 
theory of the Vlasov equation. This is the usual Boltzmann approach taking into account 
the specificities of the Vlasov equation (Casimir constraints). Therefore, the Lynden- 
Bell prediction depends on the details of the initial condition, not only on the mass and 
energy. However, violent relaxation can be incomplete: the system may not mix well 
and the collisionless relaxation may not be ergodic. In that case, other distributions can 
emerge. Among them, the Tsallis distributions (polytropes) seem to play a particular 
role in certain situations. However, they are not universal attractors. Other distributions 
can emerge as well but they are difficult to predict since they depend on the dynamics 
and on the efficiency of mixing. 

3. On longer timescales, the system evolves under the effect of "collisions", i.e. 
correlations due to graininess (finite N) effects. This gives rise to a collision term in 
the r.h.s. of the Vlasov equation. In principle, the collisional relaxation is towards 
the Boltzmann statistical equilibrium state (for stellar systems, the relaxation towards 
Boltzmann equilibrium is hampered by the problems of evaporation and gravothermal 
catastrophe). The collisional relaxation time tn(N) depends on the number of particles 
N and diverges rapidly as N — > +oo. 

4. Due to the huge timescale separation between the dynamical time tp and 
the relaxation time tn(N), the distribution function /(r, v,t) passes by a sequence of 
quasi equilibrium states which are quasi stationary states of the Vlasov equation slowly 
evolving in time due to finite N effects (graininess). In stellar dynamics, this collisional 
evolution is described by the orbit-averaged-Fokker-Planck equation. 

This scenario was put forward by Henon [231 ES] in the context of stellar dynamics 
and further discussed and improved by many authors. It has become the "standard 
scenario" of stellar dynamics [861 157] . A similar scenario was proposed by Chavanis 
[39| |4"5] for two-dimensional point vortices, by developing the analogy between the 
dynamics of stellar systems and point vortices. Inspired by these results, this scenario 
was also considered by Yamaguchi et al. [57] in relation to the HMF model. These 
authors understood many important aspects of the dynamics. However, we would 
like to complete some of their arguments in the light of more recent observations: (i) 
Yamaguchi et al. [57] argue that, during the collisional evolution, the system remains 
always close to a stable steady state of the Vlasov equation. Recently, the numerical 
study of Campa et al. [72] has demonstrated that, during the collisional evolution, 
the distribution function f(v,t) can be well-represented by a Tsallis [55] distribution 
(polytrope) with a time dependent index q(t) until it becomes Vlasov unstable and 
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triggers a dynamical phase transition from the non magnetized (homogeneous) phase to 
the magnetized (inhomogeneous) phase. Therefore, in that case, the relaxation involves 
a dynamical instability [72]. In fact, at the verge of the dynamical phase transition, 
the imaginary part of the complex pulsation Ui becomes zero and it is possible that 
the system remains marginally stable in the subsequent evolution (this would be an 
interesting point to check), (ii) Yamaguchi et al. [57] reject the possibility of Tsallis 
g-distributions [88J although they seem to play an important role in the HMF model as 
discussed in [89]. (iii) They argue that the relaxation time scales algebraically like iV 1,7 
although, as we shall see below, the picture is more complex. 

We also emphasize that the above scenario is not the only possibility and that, as 
reviewed in [90~| |3T| . the evolution can be more complex with the formation of phase 
space holes, periodic solutions, oscillations, vortex crystals etc. 

We would like now to use the kinetic theory presented in the present paper to 
interpret the results of numerical simulations that have been obtained in different 
contexts. We briefly review these numerical results, then provide new interpretations 
that can open interesting directions of research. 

Stellar systems: the scenario proposed by Henon [231 EE] has been studied by several 
authors. The theory of violent relaxation developed by Lynden-Bell [71] explains the 
isothermal core of elliptical galaxies without recourse to collisions that operate on a 
much longer timescale. However, his theory does not explain the structure of the halo 
that results from incomplete relaxation. The velocity distribution in the halo is radially 
anisotropic and the density profile decreases like r -4 (instead of the r~ 2 isothermal law). 
Theories of incomplete relaxation have been elaborated by Stiavelli & Bertin [91], Hjorth 
& Madsen [92] and Chavanis et al. [83]. Numerical simulations of violent relaxation 
have been made by van Albada [93] who compared his results with the phenomenological 
de Vaucouleur R 1 / 4 law. In the collisional regime, the evolution is first governed by 
the evaporation of high energy stars. Numerical simulations [86] show that the system 
reaches a quasi stationary state close to the Michie-King distribution that slowly evolves 
in amplitude due to the evaporation as the system loses mass and energy. The density 
follows the isothermal law r -2 in the central region (with a core of almost uniform 
density) and decreases like r _7//2 in the halo (that has anisotropic velocity distribution). 
At some point of the evolution, when the density contrast between the core and the halo 
becomes sufficiently high, the system undergoes a thermodynamical instability related 
to the Antonov [25] instability, and the gravothermal catastrophe [26] takes place. Core 
collapse has been studied numerically by Cohn [21] using the orbit-averaged-Fokker- 
Planck equation. He finds that the collapse is self-similar, that the central density 
becomes infinite in a finite time, and that the density behaves like r -2 ' 23 . When the 
system is confined in a box, instead of being free to expand, Taruya & Sakagami [91] 
observe numerically that the evolution follows a sequence of g-distributions (polytropes) 
with a time dependent index q(t) until an instability develops (that we interprete as a 
Vlasov dynamical instability in [72] )• We think that a good understanding of stellar 
dynamics has been reached by astrophysicists for a long time. 
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Point vortices: the scenario proposed by Chavanis [3H H5] , by analogy with stellar 
dynamics, has been tested numerically by Kawahara & Nakanishi [951 ESj in two very 
interesting papers. We would like to complete their discussion in the light of the kinetic 
theory presented in this paper. The kinetic theory developed in [391 HS1 SO] is valid at 
the order 0(1/N) and, consequently, the natural timescale is Nto- This is precisely 
the relaxation time that Kawahara & Nakanishi observe [96]. This seems consistent 
at first sight with the kinetic theory. However, if we consider things carefully, the 
situation appears to be more complicated. Indeed, we have seen in Sec. 13.11 that the 
collision term can vanish leading to a relaxation time larger than Ntr>. However, we have 
indicated that this cancelation occurs only for axisymmetric systems with monotonic 
profile of angular velocity. Now, in the simulations of Kawahara & Nakanishi [96J the 
system is non-axisymmetric. Indeed, in most cases, the phase of violent relaxation 
leads to a vortex crystal (see their Figs 2.b and 7). Therefore, we cannot use the 
axisymmetric kinetic equation (I68I) and we must come back to the more general kinetic 
equation (1521) . Little is known concerning this complicated equation but we expect 
that it increases the entropy and conserves energy and circulation [40J . It is not clear 
whether this equation rigorously converges towards the Boltzmann distribution but, 
since the entropy increases, we deduce that the system has the tendency to approach the 
Boltzmann distribution (it becomes more mixed). This is because there are additional 
resonances with respect to a purely axisymmetric evolution. Since the proper timescale 
of this kinetic equation is Ntu, the kinetic theory of non-axisymmetric flows may 
explain the numerical observations of Kawahara & Nakanishi [96]. There remains, 
however, a difficulty. Indeed, in their Fig. 2a, violent relaxation leads to an apparently 
axisymmetric distribution with a presumably monotonically decreasing profile of angular 
velocity. In that case, there should be no resonance and the relaxation time should be 
larger than Nto contrary to what is observed numerically. We speculate that, in this 
case, the system is not really axisymmetric but that there exists "structures" in the 
flow (analogous to Dawson's vortices in neutral plasmas under a strong magnetic field 
[97J). This claim is corroborated by the fact that Kawahara & Nakanishi [96] observe 
anomalous diffusion which means that the evolution of the flow is more complex than 
it seems. These results have been confirmed recently by Yoshida [9H] who showed that 
the mean square displacements exhibit a universal time dependence ([r(t) — r(0)] 2 ) ~ t M 
with /i = 1.75 ± 0.1. These authors attribute anomalous diffusion to occasional long 
jumps of the particles convected by long living large vortices. If there exists such 
"structures", they may induce new resonances and explain the relaxation towards 
Boltzmann equilibrium on a timescale Nto- This may also imply that there is no 
contradiction with the kinetic theory concerning the diffusion process. The diffusion 
coefficient (!79|) and the Fokker-Planck equation (|82|) derived in [39l H5] are only valid for 
axisymmetric flows without "structures" . In more general cases, the diffusion equation is 
given by Eq. (142) of [ID] which involves derealization in space and time. It is plausible 
that these non ideal effects lead to anomalous diffusion for short times. Normal diffusion 
should be recovered for larger times (anomalous diffusion could also be a finite size effect 
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and it could disappear for large N) . As indicated by Kawahara & Nakanishi [9B] , this is 
precisely what happens for neutral plasmas under a strong magnetic field in relation to 
the Taylor-McNamara [99J diffusion. Of course, this scenario demands to be confirmed 
but it seems plausible enough to reconcile theory and observations. More work remains 
to be done in order to better understand these results. In this respect, we may note 
the recent study of Sano |100j on the kinetic theory of point vortices. He confirmed 
the results of Chavanis [391 SO] (leading to Landau-type kinetic equations) and derived 
a kinetic equation corresponding to the Lenard-Balescu equation for plasmas using the 
theory of the Fredholm integral equation. We may also note that the exponential tail of 
the velocity distribution observed by Kawahara & Nakanishi [96] is consistent with the 
results reported by Chavanis & Sire |101j when the vortices have a finite core, which 
is equivalent to putting a restriction on the point vortex distance in the initial state. 
Therefore, this observation (exponential tail) is not inconsistent with kinetic theory. 

HMF model: the scenario proposed by Yamaguchi et al. [57] has been checked 
numerically by these authors and by Campa et al. [581 172] . The relaxation time has 
been discussed by Bouchet & Dauxois [52] and Chavanis [31] in the light of the kinetic 
theory. These authors argue that, for spatially homogeneous systems, the relaxation 
time is larger than Ntp since there is no resonance. Bouchet & Dauxois [52] claim that 
this result is consistent with the timescale N lJ numerically observed by Yamaguchi et 
al. [57]. Chavanis [31] claims that this is consistent with the timescale e N observed by 
Campa et al. [58]. We would like to reconsider and complete their discussions. In the 
simulation of Yamaguchi et al. [57], that is performed for subcritical energies U < U c , 
the system does not remain spatially homogeneous so that the homogeneous kinetic 
theory cannot be advocated until the end. Alternatively, Campa et al. [58] perform 
supercritical simulations U > U c , in which the system remains permanently spatially 
homogeneous, and find that the relaxation time scales like e N . Therefore, when the 
system is homogeneous the relaxation time is exponentially long! This is in agreement 
with the kinetic theory but this shows that the relaxation is not due to just three-body, 
four-body, n-body... correlations (as we could naively believe), that would imply a 
relaxation time scaling as N 2 , N 3 , iV n_1 ... All the correlations functions seem to matter! 
Therefore, the kinetic theory of spatially homogeneous systems cannot apparently be 
truncated at some power of 1/N. Using the analogy between the point vortex model 
and the HMF model, we may conjecture that the relaxation time for an axisymmetric 
distribution of point vortices with monotonic profile of angular velocity scales like e N 
(although the possibility that such distributions never achieve statistical equilibrium 
cannot be rejected). When the system is spatially inhomogeneous, the relaxation time 
is reduced with respect to e . This leads to a reinterpretation of the numerical results. 
Yamaguchi et al. [57] performed two types of simulations: (i) the first one, corresponding 
to unstable stationary Vlasov states, leads to a spatially inhomogeneous evolution. In 
that case, the collisional relaxation times scales like N which is the proper scaling 
of the kinetic theory for spatially inhomogeneous systems. Indeed, there exists many 
resonances that "push" the system towards the Boltzmann distribution (see Sec. H]). The 
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question whether the system really reaches the Boltzmann distribution, or a distribution 
close to it, remains open, (ii) The second, which starts from Vlasov stable stationary 
states, leads to spatially homogeneous structures. In that case, there is no resonance and 
we should expect a relaxation time scaling as e . However, Campa et al. [72] argue that 
the relaxation towards Boltzmann equilibrium involves a dynamical phase transition 
that is triggered when the homogeneous phase becomes Vlasov unstable. In that case, 
the system becomes spatially inhomogeneous and resonances appear. This considerably 
accelerate the relaxation with respect to the exponential scaling and leads to a non 
trivial (but probably non universal) scaling N 1,7 . Therefore, in that case, the timescale 
cannot be solely understood in terms of the kinetic theory of homogeneous systems, 
although this theory explains why the scaling is larger than N . Morita & Kaneko 
[102] . considering an initial condition with U < U c and Mq = 1, found a relaxation 
time tf> ~ Nto- In their simulations, the system is always spatially inhomogeneous 
(the magnetization has an oscillatory behavior) and the relaxation time is consistent 
with the natural scaling of the kinetic theory for inhomogeneous distributions. Their 
results may be consistent with the general kinetic equations ([2]) and fllOip although this 
is of course difficult to check. Let us finally consider the early simulation of Latora 
et al. |103j which starts from an out-of-equilibrium initial state with magnetization 
Mq = 1. Violent relaxation leads the system to a spatially homogeneous QSS which 
is followed by a collisional relaxation towards the Boltzmann distribution. They find 
a relaxation time scaling like Ntjj. This timescale seems to be inconsistent with the 
kinetic theory of spatially homogeneous systems. However, as in the case of 2D point 
vortices, we speculate that there exists "phase-space structures" that make the system 
spatially inhomogeneous and induce additional resonances. These inhomogeneities may 
reduce the timescale of relaxation to the natural timescale Ntp. These "phase-space 
structures" are indeed present in the simulations of Latora et al. [103] and they seem 
to be responsible for anomalous diffusion in a way similar to what happens in 2D point 
vortex dynamics [96], [98] . Again, this scenario demands to be confirmed but it seems 
plausible enough to reconcile theory and observations. We must, however, be very 
careful because these striking features (phase-space structures and anomalous diffusion) 
may be due to finite size effects [104[ 1105] and disappear for N — > +oo. Also, diffusion 
may be anomalous for short times and normal for large times. This may be a way to 
reconcile the approaches of Rapisarda & Pluchino [106] and Bouchet & Dauxois [52] who 
studied temporal correlation functions and anomalous diffusion. Indeed, their studies are 
based on very different arguments: Bouchet & Dauxois [52] assume that the background 
distribution is spatially homogeneous and use standard kinetic theory based on Fokker- 
Planck equations. Alternatively, Rapisarda & Pluchino [106] assume that anomalous 
diffusion is due to the presence of phase-space structures and use nonlinear Fokker- 
Planck equations. We believe that these two approaches are not antagonistic but that 
they apply to different regimes or, maybe, different numbers of particles. 

In conclusion, we have developed a kinetic theory of systems with long-range 
interactions based on the generalized (and not well-known) kinetic equations ([2]), (1521) 
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and (EH}. A specificity of our approach is to remain in physical space and take 
into account spatial inhomogeneity which is the norm for systems with long-range 
interactions. This is the main difference with kinetic theories of systems with short- 
range interactions and neutral plasmas that are spatially homogeneous. Interestingly, 
the homogeneous kinetic equation can be easily derived from the generalized kinetic 
equation with simple calculations. This provides a new derivation of Landau-type 
equations and also offers a unified description of the kinetic theory of stellar systems, 
2D point vortices and the HMF model. As an application, we have argued that 
the relaxation time of spatially inhomogeneous ID systems (or non axisymmetric 
2D flows) should scale as Nto which is the natural timescale of the kinetic theory. 
Indeed, spatially inhomogeneous systems involve many resonances that allow the 
relaxation (or the approach) of the system towards Boltzmann equilibrium. By contrast, 
for one dimensional homogeneous systems (or axisymmetric flows with monotonic 
profile of angular velocity), the resonances disappear and the relaxation timescale is 
considerably increased. We have speculated that, when the system remains permanently 
homogeneous, the natural timescale is e N t£>. In the HMF model, this corresponds 
to supercritical energies considered by Campa et al. [58]. If the QSS resulting 
from violent relaxation is homogeneous but, later evolving under the development of 
correlations, undergoes at some point a dynamical phase transition (Vlasov instability) 
from the homogeneous phase to the inhomogeneous phase, the relaxation time should be 
intermediate between Ntp and e N tp. For example, it can scale algebraically as N 5 t^. 
This corresponds to subcritical energies considered by Yamaguchi et al. [57] . They 
find 5 = 1.7 but this scaling may not be universal. We have motivated our arguments 
by a detailed comparison with existing numerical simulations and found interesting 
agreements. This shows that kinetic theory can be a good approach to understand the 
dynamics and thermodynamics of systems with long-range interactions but that many 
works remain to be done to obtain a complete description of spatially inhomogeneous 
systems. 

Appendix A. Calculation of K^ u 

In d — 3, we introduce a spherical system of coordinates where the z axis is taken in 
the direction of w. Then 

K^ U = 7r(27r) 3 m 2 / k 2 dk / sin9d9 d<p k^k u 5(kw cos 9)u{k) 2 . (A.l) 
Jo Jo Jo 

Using k x = k sin 9 cos 0, k y = k sin 9 sin <fi and k z = kcos9, it is easy to see that only 

K xx , K yy and K zz can be non- zero. The other components of the matrix vanish 

by symmetry. Furthermore 

K xx = K yy = vr 2 (27r) 3 m 2 / k 2 dk sin 9 d9 k 2 S(kw cos 9)u{k) 2 sin 2 9. 

Jo Jo 

(A.2) 
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Using the identity 5(Xx) = ■rhS(x), we get 

1 r+oo rn 

K xx = K yy = 7r 2 (27r) 3 m 2 - / k 3 u(k) 2 dk sin 3 0<J(cos0) dO. (A.3) 

w Jo Jo 



w Jo 

With the change of variables t = cos 6, we obtain 



I r+oo r + l 

K xx = K yy = vr 2 (2vr) 3 m 2 - / k 3 u(k) 2 dk (1 - t 2 )5(t) dt, (A.4) 

w Jo J-i 

so that finally 

1 r+oo 

K xx = K yv = 8vr 5 m 2 - / k 3 u(k) 2 dk. (A.5) 
w Jo 

On the other hand 

K zz = 27r 2 (27r) 3 m 2 - / k 3 u(k) 2 dk / sin 9 cos 2 95(cos 9) d9 = 0. (A.6) 
w Jo Jo 



In conclusion, we obtain 



K» v = — [ 5^ - ) , (A.7) 



w \ w 

with 

r+oo 



r+oo 

K = 8n 5 m 2 k 3 u(k) 2 dk. (A.8) 
Jo 

In d = 2, we introduce a polar system of coordinates where the x axis is taken in 
the direction of w. Then 

r+oo r2ir 

= n(2ir) 2 m 2 / kdk / d9 k»k u 5 {kw cos 9)u{k) 2 . (A.9) 
Jo jo 

Using k x = k cos 9, k y = k sin 9, it is easy to see that only K xx and K yy can be non-zero. 
The other components of the matrix K^ u vanish by symmetry. Furthermore 

r+oo r2ir 

K yy = n(2ir) 2 m 2 kdk d9 k 2 5(kw cos 9)u{k) 2 sin 2 9, (A.10) 
jo jo 

or, equivalently, 

1 r+oo /*7T 

K yy = 2n(2n) 2 m 2 - / k 2 dku(k) 2 d9 5(cos 9) sin 2 0. (A.ll) 
Jo jo 

With the change of variables t = cos 9, we obtain 

1 r+oo r+l 

K yy = 27r(27r) 2 m 2 - / k 2 u(k) 2 dk y/T^¥S(t)dt, (A.12) 
w Jo J-l 

so that finally 



1 r+oo 

K yy = 8n 3 m 2 - k 2 u(k) 2 dk. (A.13) 

U7 JO 

On the other hand 

1 r+oo r2n 

K xx = n(2ir) 2 m 2 - k 2 dku(k) 2 d9S(cos9) cos 2 9 = 0. (A.14) 
w Jo Jo 



In conclusion, we obtain 



K^ = -^-^f), (A.15) 



w \ w 

with 

r+oo 



r+oo 

K = 8Tr 3 m 2 k 2 u{k) 2 dk. (A.16) 

JO 
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Appendix B. Calculation of x( r ; r i) 



For the ordinary potential u(\r — rx|) = — j- In |r — ri|, we have 
u(r, r±, 4>) = — — ln(r 2 + r\ — 2rr\ cos</>), 

47T 

where <fi = 9 — Q\. Its Fourier transform on the angles can be written 

1 f*, . 5m0 



(b.i; 



u m {r, r lj 



471-2 



/"7T m Q 

/ ln(l — A cos 0) cos(m0) dcf) ^-ln(r 2 + r 2 ), (B.2) 

JO 4:71 



with A 



^2+l-i ■ For m ^ 0, we have the identity 



In (1 — A cos <b) cosfmc 



7T / 1 



m \ A 



Therefore 



u m (r, n 



^2 _|_ ^,2 |t"^ 7"^ | \ 



1 />< 
47r|m| 1 r- 



\m\ 



(B.3) 



(B-4) 



47r|m| \ 2rr\ J <±n\ni>\ \/> / 

where r< (resp. r>) is the min (resp. max) or r and r±. For m = 0, we have the identity 



r In (1 - A cos <p) d<p = tt ln(l + Vl - A 2 ) - tt In 2. 

JO 



Therefore 



A . . 1 , ( r 2 + r? + |r 2 — rf I 
«o(r, ri) = -— In I 



In r> . 

2tt > 



Combining the previous results, the potential of interaction can be written 

m\ 



u (r,n,0) = --^lnr> + — ^ — ( — ] 
2vr 47r^ |m| V?W 



rn^O 

The function defined by Eq. ( |69l) takes the form 

2m 

x(r, n) 



1 g 1 lr< 



m=l 









In 


1 - 


(?) 









In particular 



1 +0 ° 1 



lnA, 



m=l 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



where In A is a Coulombian logarithm that has to be regularized appropriately (see, e.g. 
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